Enhanced Orbital Degeneracy in Momentum Space for LaOFeAs 



Hai-Jun Zhang, Gang Xu, Xi Dai and Zhong Fang 
Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China; 

(Dated: April 17, 2008) 

The Fermi surfaces (FS) of LaOFeAs (in fc z =0 plane) consist of two hole- type circles around V 
point, which do not touch each other, and two electron-type co-centered ellipses around M point, 
which are degenerate along the M-X line. By first-principles calculations, here we show that addi- 
tional degeneracy exists for the two electron- type FS, and the crucial role of F-doping and pressure 
is to enhance this orbital degeneracy. It is suggested that the inter-orbital fluctuation is the key 
point to understand the unconventional superconductivity in these materials. 

PACS numbers: 74.7D.-b, 74.25. Jb, 74.25.Ha, 71.20.-b 



The discovery of superconductivity in LaOFeAs fam- 
ily of compounds is really exciting. Very soon 
after the discovery of up to 26 K superconduc- 
tivity in La(Oi_ :c F a: )FeAs Q, the T c =43 K in 
Sm(Oi_ x F x )FeAs 0, T C =4T K in Ce(Oi_ a; F :c )FeAs 
and T c =52 K in PrOi_ a: F a; )FeAs were reported. More 
and more experimental and theoretical studies suggest 
that the superconductivity found here is unconventional 
non-BCS type 0,0,0]. Several possible pairing mech- 
anisms have been discussed theoretically 0, [t| , however 
we show here, based on first-principles calculations, that 
the key is the inter-orbital fluctuation due to orbital de- 
generacy at the two electron-type Fermi surfaces, which 
is much enhanced by either F-doping or pressure. 

BE, 



The electronic structure 



12, 131 of LaOFeAs is 



quasi two dimensional and very similar to typical semi- 
metal. Namely there are three hole-type Fermi sur- 
faces (FS) around the T-Z line, and two electron-type 
FS around the M-A line of the Brillouin zone (BZ). The 
existence of five FS makes the problem complicated, how- 
ever the electronic structure around Fermi level can be 
qualitatively understood in a simple way using three or- 
bitals per Fe site (d xy , d zx , d yz ). Here we define the local 
coordinates x and y to be 45° rotated from the a and b 
directions of the crystal lattice respectively. 

Let's focus on the k z =0 plane. There are two co- 
centered circles around the T point, which are hole-type 
FS and mostly come from the d yz and d zx states of Fe. 
In addition, there are two co-centered ellipses around the 
M point, which are electron- type FS and mostly come 
from the d xy and d yz / zx states of Fe (see Fig.l). The 
anisotropy (or the shape) of the ellipses is related to the 
hybridization between d xy and d yz (or d zx ). The stronger 
the hybridization, the more distortion of the ellipses. The 
projected orbital characters of the electron-type FS are 
shown in Fig.l. One important fact to be noted here is 
that the two hole-type circles around T do not touch each 
other, however, the two electron-type ellipses around M 
are degenerate along the X-M line (see Fig.l), and this 
degeneracy is protected by the crystal symmetry of the 
system. However, we will show here that additional dc- 
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FIG. 1: Upper panels: Schematic plots for the Fermi surfaces 
of LaOFeAs. The left panel is for pure LaOFeAs, while the 
right panel is for the case with about 20% F-doping. The 
large square (solid lines) is the BZ of one Fe cell, while the 
small diamond (dashed lines) is the BZ of two Fe cell. Please 
note the interchange of the orbital characters of two electron- 
type FS after F-doping. Lower panel: The projected orbital 
characters of states at FS, which is defined as |(a|\Pfc)|. Here 
a are five local Fe-3d orbitals, Vl/j. are Bloch states around one 
of the ellipses, and angle 6 is defined in the inset. 



generacy exists for the two branches of electron- type FS, 
and the crucial roles of F-doping or pressure is to enhance 
this degeneracy. 

As discussed recently by several papers HQ, signifi- 
cant FS nesting exists between the hole and the electron 
type FS for pure LaOFeAs. After shifting the hole FS 
around T by a vector q—(n, 7r,0), it will largely over- 
lap with the electron FS around M. This leads to spin- 
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FIG. 2: Calculated Fermi surfaces of La(Oi_ a; F a ;)FeAs around 
the M point for different F-doping levels. Only part of the 
BZ is show and the hole-type Fermi surfaces around V point 
is not included. From top to bottom are, (a) no F-doping; (b) 
5% F-doping; (c) 10% F-doping; (d) 15% F-doping. The two 
ellipses are almost degenerate for 10% F-doping. 



density-wave (SDW) instability, and is the origin of the 
150 K anomaly observed by transport, susceptibility, and 
optical measurements [l4j . It is now realized that the 
SDW instability actually competes with superconductiv- 
ity, and super-conducting state can be realized only after 
the SDW transition is suppressed 0, [3| . Electron dop- 
ing by F is an efficient way to suppress the nesting effect 
(the SDW instability) , because the shifting of Fermi level 
enhances the mismatch between electron and hole FS (as 
shown in Fig.l). It is also true that the superconductiv- 
ity was observed in LaONiAs 15j, where nesting effect 
discussed above is irrelevant. All those information sug- 
gest that the nesting between the electron and hole FS 
can not be the key issue for the superconductivity. 

Any role else played by F-doping in achieving high T c 
superconductivity? This is the main point we would like 
to address here, based on quantitative first-principles cal- 
culations. We use plane wave pseudo-potential method 
and the generalized gradient approximation (GGA) for 
exchange-correlation potential, the same approach as we 
reported before As discussed above, two hole- type 
FS are co-centered circles around T point. They have 
different size and do not touch each other. However, 
two electron- type FS are ellipses around M point. They 
have the same size, but are rotated by 90° from each 
other. As the result, they overlap (or degenerate) at spe- 
cial k-points along the X-M line of the BZ (see Fig.l 
for schematic understanding). The size of the electron- 
type FS will increase with F-doping, however, most in- 
terestingly the shape of two ellipses is also significantly 
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FIG. 3: Calculated band structure of LaOFeAs for (a) zero 
pressure; (b) 30 GPa pressure; (c) fixed volume but reducing 
the c/a ratio by 5%. The crossing of two eletron-type bands 
is indicated by red circle. 



changed by F-doping, as shown in Fig. 2. With the in- 
crement of F-doping, the anisotropy of ellipses is sup- 
pressed (i.e, they are more and more circle-like), as the 
result the orbital degeneracy right at the two electron- 
type FS is enhanced. By about 10% F-doping, the two 
ellipses become very close to circles and almost exactly 
overlap. However, with further F-doping, two ellipses are 
distorted again. 10%-F is about the optimal doping level 
found experimentally for this family of compounds, and 
our calculations show that at this doping level the orbital 
degeneracy in the momentum space is the highest. 

More careful analysis about the band structure near 
the FS show that additional degeneracy exists at the two 
electron-type FS, which can be tuned by F-doping. This 
degeneracy can be understood from the calculated band 
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structure as shown in Fig. 3(a). Looking at the M~r di- 
rect, there are two dispersive electron- type bands cross- 
ing the Fermi level. However, the two bands disperses 
at different velocity and they cross each other at about 
85 meV above the Fermi level. It is this band-crossing 
which leads to the enhancement of FS degeneracy dis- 
cussed above. This band-crossing is accidental, but it 
is guaranteed by the different dispersion velocities of the 
two bands, and it is a common feature as seen in this fam- 
ily of compounds [T3| ■ One of the bands has mainly char- 
acter of two-dimensional d xy orbital, which is strongly 
dispersive in the plane; another one of the bands is the 
linear combination of one-dimensional orbitals d yz and 
d zx orbitals. The key point is that the low-lying d xy 
band (at M point) disperses faster than the high-lying 
dyz/zx band. By up-shift of the Fermi level (F-doping), 
the two electron-type FS will become more and more de- 
generate, however beyond 85 meV above the Fermi level, 
the two FS interchange the orbital character and start to 
deviate from each other (as shown in Fig.l and Fig. 2). 

In the following we discuss the effect of pressure and 
point out that pressure can play a similar role as F- 
doping. For this purpose, we fully optimize the crystal 
structure (including internal coordinates) using the ex- 
perimental space group PA/nmm for each fixed volume. 
Then the electronic structures are calculated using the 
optimized structure. Fig. 4 shows the calculated total en- 
ergies, c/a ratios, and internal coordinates as function of 
volume. The general tendency of lattice behavior under 
pressure is quite smooth and the calculated bulk mod- 
ulus is about 90 Gpa, suggesting that the compound is 
relatively easy to be compressed. Other details about 
the effect of pressure will be discussed elsewhere, here we 
concentrate on the electronic structure. Fig. 3(b) shows 
the calculated band structure of LaOFeAs under pressure 
(about 30 Gpa). It is clearly seen that the crossing be- 
tween the two electron-bands (along M-r line) is shifted 
downward to be less than 1 meV above Fermi level. In 
the meanwhile, the size of the electron-type FS is en- 
larged. All these indicate the same tendency as discussed 
above for F-doping, namely the orbital degeneracy of two 
electron-type FS is enhanced by pressure. In addition to 
the uniform pressure, we also found that this degeneracy 
can be enhanced by uniaxial pressure along the c axis. As 
shown in Fig. 3(c), the band crossing is shifted downward 
to be about 7meV above Fermi level by 5% shorting of 
c/a ratio with fixed volume. 

The high orbital degeneracy at the FS will lead to 
strong inter-orbital fluctuation and may lead to inter- 
orbital pairing. In such case, it will be very easy to un- 
derstand why the T c is enhanced initially by F-doping [lj 
and pressure and is suppressed beyond the optimal 
doping level about 10%. Based on this inter-band pair- 
ing picture, the possible pai ring symmetry has been dis- 
cussed in our parallel paper ljfl . Here we will show that 
the experimentally observed T c versus F-doping phase di- 
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FIG. 4: The calculated total energy and optimized structure 
parameters as function of volume. Here the internal structural 
parameter d z are dinned as the inter-layer distance along c 
axis. 



agram can be explained semi-quantitatively by our sce- 
nario. For this purpose, we construct a minimum two 
bands effective Hamiltonian to describe the electron-type 
FS. 



As we have discussed in our parallel paper [17|, the 
strong inter-orbital ferromagnetic fluctuation will induce 
an even parity, orbital singlet and spin triplet paring 
state. This tendency will be strongly enhanced by the 
large orbital degeneracy at the Fermi surface. Since the 
electron pairing happens in momentum space, the main 
contribution to the paring comes from the area near the 
FS. Therefore for the purpose of illustration, we can use a 
very simple two band model on the square lattice, which 
not only gives the right shape but also the density of 
states near the FS for different F doping. The Hamilto- 
nian reads, 



k'm 



kk J 



where a = 1, 2 are orbital indicies, and 



(1) 



(2) 



with 7^ = cosfc x + cosfcy,7^. = cos(k x + k v ) and 7^ = 
cos(k x - k y ). 

By tuning the hopping integrals t, t\ and t%, we can ob- 
tain very similar FS with the calculated results from first- 
principles for different doping. We then solve the above 
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sure is to enhance the orbital degeneracy, which leads 
to strong inter-band fluctuation. Combining with exper- 
imental observations, we suggest that the inter-orbital 
pairing is the key to be considered for the mechanism of 
superconductivity in this compound. The same orbital 
degeneracy is also present in LaONiAs, and this will be 
discussed in our separated paper (lij . 
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FIG. 5: The super-conducting gap as the function of doping 
concentration for LaOi-^F^FeAs. The positive x is defined 
for electron doping, while negative x for hole-doping. The 
inter-band ferromagnetic coupling J is chosen to be 0.45 eV. 



Hamiltonian by mean field decoupling, 
detail in our previous paper 
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as described in 
The spacial pairing 
symmetry is chosen to be s-wave and k dependence of 
is also ignored for simplicity. For each doping, we first fit 
the FS to the GGA results by tuning the above hopping 
integrals and then calculate the super-conducting gap by 
mean field approach. The results are plotted in Fig. 5. Al- 
though our calculation is very rough, it still indicate two 
important features of our theory. First, the superconduc- 
tivity in the hole doped side and electron doped side is 
from the same origin. As pointed out previously [14| , for 
the pure compound (i.e. without doping), there is an ex- 
tra SDW instability due to the nesting effect between hole 
and electron FS, which kills the superconductivity. Both 
electron and hole doping will remove the FS nesting effect 
and superconductivity appears. Second the strongest su- 
perconductivity tendency appears on the electron doping 
side, with the doping concentration around 10%. Both 
of these predictions seem quite consistent with the cur- 
rent experimental data. On the other hand, care must 
be taken if both F-doping and pressure are applied. Our 
prediction is that T c will be enhanced by pressure for 
samples below optimal F-doping (about 10%), however 
it is opposite for overdoped samples. 

In summary, by first-principle calculation, we show 
that the crucial role played by either F-doping or pres- 
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